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Abstract 

An exact solution is given for a two-dimensional model of a Coulomb gas, more general 
than the previously solved ones. The system is made of a uniformly charged background, 
positive particles, and negative particles, on the surface of a sphere. At the special value 
r = 2 of the reduced inverse temperature, the classical equilibrium statistical mechanics is 
worked out : the correlations and the grand potential are calculated. The thermodynamic 
limit is taken, and as it is approached the grand potential exhibits a finite-size correction 
of the expected universal form. 



KEYWORDS : Coulomb gas ; solvable models ; finite-size corrections. 



LPTHE Orsay 95-48 
June 1995 



1 Department of Mathematics, University of Melbourne, Parkville, Victoria 3052, 
Australia. E-mail : matpjf@maths.mu.oz.au 

2 Permanent address : Laboratoire de Physique Theorique et Hautes Energies, Univer- 
site de Paris XI, batiment 211, 91405 Orsay Cedex, France (Laboratoire associe au Centre 
National de la Recherche Scientifique - URA D0063). E-mail : janco@stat.th.u-psud.fr 



2 



1. INTRODUCTION 

At the special coupling V = q 2 /ksT = 2, where q denotes the magnitude of the 
charges, T the temperature, and ks Boltzmann's constant, the classical two-dimensional 
one-component and two-component plasma systems of point particles can be solved 
exactly in a number of different geometries. In the one-component case^ 1-3 -* the calcula- 
tion is based on the Vandermonde determinant identity and is performed in the canonical 
ensemble, while in the two-component case^ 4-7 ) it is based on the Cauchy determinant 
identity and is performed in the grand canonical ensemble. Also, in the one-component 
case the classical Boltzmann factor at V = 2 is isomorphic^ to a squared wave function 
for N non-interacting fermions in two-dimensions, subject to a uniform magnetic field and 
confined to the lowest Landau level. The grand partition function of the two-component 
plasma at V = 2 is isomorphic^ 5 ) to the generating functional of the two-dimensional free 
Fermi field. 

In Sections 2 and 3 of this paper we unify the exact calculation of the particle cor- 
relations for one-component ( 2 ) and two-component^ plasma systems on a sphere. This 
is achieved by introducing a new formalism based on a determinant identity^ 9 -* which 
interpolates between the Vandermonde and Cauchy identities. In fact the formalism is 
of general applicability in the sense that it is not restricted to the sphere domain. The 
exact calculation is performed in an ensemble with a uniform background charge density of 
given total charge — Nq, a given number N of positive particles, and in addition a variable 
number M of pairs of positive and negative particles which is controlled by a fugacity £• 
In the limit £ — > the canonical ensemble of the one-component plasma is reclaimed, while 
in the limit N — > it is the grand canonical ensemble of the two-component plasma which 
is reclaimed. 

In the infinite plane system the correlations for the two-component plasma plus back- 
ground system at V = 2 are already known from the work of Cornu and Jancovici/ 10 -' who 
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used the formalism for the two-component plasma in a field^ 5 -' to perform the calculation. 
This suggests an alternative approach, considered in Section 4, to the system on a sphere : 
at the north pole fix a point charge Nq and then use the formalism of ref. 5. Of course the 
north pole becomes highly singular, as we expect that N particles of negative charge will 
accumulate about this point. However, due to the screening of the charge-charge correla- 
tions, we would expect this not to affect the correlations between charges away from the 
north pole. In fact, we find that this method reproduces the exact correlations of Section 3. 

The grand potential Q# can be calculated by integration with respect to ( of the 
expression for the density of negative charges. The large- R asymptotic expansion of is 
of particular interest as for a Coulomb system on a sphere it has been predicted that^ 6 ) 

fi fl -f2oo~^lnii , (1.1) 

independent of the detail of the plasma. In Section 5, (1.1) is verified. 

In the Appendix, it is shown how the partition function for the one-component plasma 
on a sphere can be calculated within the formalism introduced in Section 2. 

2. FORMALISM BASED ON A DETERMINANT IDENTITY WHICH 
INTERPOLATES BETWEEN THE VANDERMONDE AND CAUCHY 
IDENTITIES 

In this Section the generalized (with position-dependent fugacities) grand canonical 
partition function for the two-component plasma plus background system at V = 2 will be 
expressed as the determinant of an integral operator. This allows the one-component and 
two-component plasma systems to be solved using the same technique. We will assume 
that the domain is the surface of a sphere for calculating the particle correlations in Section 

3. However, the formalism could also be carried through in the other domains for which 
the one- and two-component plasmas are solvable, in particular the disk and semiperiodic 
boundary conditions. 



For two charges q and q' on the surface of a sphere of radius R, separated by an 
angular distance ip, the pair potential is^ 2 ) 

<j)(^) = -qq'ln[2Rsm(tjj/2)/£} , (2.1) 

where £ is an arbitrary length scale (to be taken as unity) and 2Rsm(ip/2) is equal to the 
chord length between the particles. 

Representing the particle positions by spherical polar coordinates 9 and <p (9 is the 
angle from the north pole and cp is usual other polar angle), and similarly 9' and <p', we 
have the identity*- 2 -* 

(2.2a) 



sin(^/2) = \a(3' - a'(3\ = \(3f3' 
where 



a a 



a = cos -e^ /2 , = -i sin - e "^/ 2 (and similarly a', (3') (2.2b) 

The coordinates a and f3 are referred to as the Cayley-Klein parameters. 

Consider now a system of N + M particles of charge q with coordinates specified 
by Cayley-Klein parameters a'-, /?'• (j = 1, • • • , N + M), M particles of charge — q with 
coordinates specified by Cayley-Klein parameters otj, (3j (j = 1,---,M), and a uniform 
neutralizing background. From (2.1) and (2.2), the full Boltzmann factor (including the 
constants from the background-particle and background-background interactions) is read- 
ily calculated as 

1 \ r(7V/2+M) M N+M 

" ' /A H'M- iv+N) n 



/ 1 \ r(7V/2+M) M N+M 

W(N + M,N;T)=(J R j e rN2/4 JJ |^| _r(1+iV) l^f^- 1 ) |D| r (2.3a) 



where 

n («fc-«j) n ( u fc- v j) 

l<j<fc<M l<j<fc<7V+M , . 

D := WnTm ( 2 - 36 ) 

n n ~ 

j=l k=l 



with 



<± and V ..-?L 
and c, . s . 



(2.3c) 



Thus the generalized grand canonical partition function, with one-body potentials with 
Boltzmann factors a(6>', <p') and b(6, (p) coupling to the positive and negative charges re- 
spectively, and with N fixed and M summed over, is given by 



i c 



2M 



Z r (a,b) = A Nr ^ m {N + M )\ \2R 



1 



TM 



R 



4M 



'N+M k 

x I n / d^sin^[sin(^/2)] 



2tt 



dtp't a{0' e , tp' e ) 



x 



M «7r »2tt \ 

IT / d^sin^[sin(^/2)]- r(1+JV) / d<p e b(9 e ,<p e ) ) \D\ r (2.4a) 



where 



/ 1 \^ r / 2 



(2.46) 



In the case N = (two-component plasma) D can be written as a Cauchy determinant, 
while in the case M = (one-component plasma) D can be written as a Vandermonde 
determinant. In the general case we have^ 9 ) 



D = (_l)M(M-l)/2 detA 



(2.5a) 



where 



A 



1 1 

Ul Um+JV 

^1 V M+N 
1 1 

Wl— Vi til- V M +N 
1 1 



(2.56) 



and consequently 

where O denotes the (M + N) x (M + N) zero matrix. 

Consider now the generalized grand partition function (2.4a) for T = 2 with the 
substitution (2.6). Suppose we discretize the integrals in (2.4a) by making the replacements 

K 1 



P.-K -i -TV 1 / 

\ dOtw*.-)^— y, / 



7i(m e - vx) 
Tv(m' e - i/j) 



27T i K 



K 2 



where < i^i, v[, v 2 , v' 2 < 1. The domain available to the particles now consists of separate 
sublattices for the positive and negative charges {v\ , • • • , v 2 are to be chosen so that no 
two lattice points overlap). We impose the constraint that the number of lattice points 
available to the positive charges, K[K 2 , is N greater than the number available to the 
negative charges, K\K 2 . The replacements (2-7) are exact in the limit Ki, ■ ■ ■ , K' 2 — > oo 
so we can write 
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K-lK 2 



K lt Kf-+oo K 2 ,K^oo ^ o M\(M + N)\ V 2 



E 2 [a,b] = A N2 lim lim 



2M 



x II \jc E h^-MIanirK-^r^ 

£=1 [ 1 m^=l 



2 »;=i 



X 



1 * 2 1 

x — ^ 6(7r(m^ - vij/K^ 2n(n t - v 2 )/K 2 ) \ 

2 ne=l ) 



X det 



O iA 
iA+ O 



(2.8) 



where the coordinates implicit in A are taken as lattice points according to the prescription 
(2-7). 

We now make the crucial observation that expression after the limits in (2.8) is pre- 
cisely the expanded form of a single determinant : 



S 2 |a, b] = lim lim det ( t' + 



iC O 



(2.9) 



where 1' is a IK'^K'^ x ^K[K 2 diagonal matrix with the first N diagonal entries zero and 
the rest one, O is the K[K 2 x K[K 2 zero matrix and 



i p— l 



(Cfl/2)b,- fc 

K 'l K '2( U jk-Vj/ k /) 



(jk) = (ll),---,(K 1 K 2 ) 

0'fc')=(ii),-.-,(if;^) 



(2.10a) 



C 



J j>k' 



Clj'k' 



0'fe') = 
p=l,." 



(11),- 



{CR/2)a rk , 
K 1 K 2 (u jk - Vj'k') 



(/*')= 

(jk) = (ll),---,(K 1 K 2 ) 



(2.106) 
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(the upper subscript on the matrix symbols [ ] labels the rows while the lower subscript 
labels the columns), where 

b jk =sin7r(j - v{]/K\ (shi7r(j - u 1 )/2K l )~ 2( - 1+N) 

x b (7r(j - vJ/Ku 27r(k - u 2 )/K 2 ) (2.11a) 



a rv =sm7r(j' - u[)/K[ (smvtf - v[) /2K[y 2{1 - N) 



(2.116) 



and Ujk denotes the coordinate u = a/ (3 with 9 and cp therein evaluated at the lattice point 
9 = w(j — vi)/Ki, ip = 2ix(k — V2)jKi (and similarly for the meaning of Vj'k> in terms of 
v = a' 1(3'). 

From (2.10) and (2.11) we see that the limits in (2.9) can be taken to give 



E 2 [a,b] = A N2 det(X) 



(2.12) 



where the operator X acts on vectors 



and is defined by 



where 



V>:= 



ai 



m<p) 



A N 

F(6,<p) 
G(6,<p) 



Ai := i 



(2.13) 



(2.14) 



p7T p2lY 

/ d9sm9 dip[sm(9/2)} N ~ 1 [cot(9/2)e lv ] J ~ 1 g(9,ip) (2.15a) 
Jo Jo 
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F(e,<p) =m<p) - ^ r d9 1 sm9 1 [sin^/2)]^- 1 [sin^)]" 1 "^ 

z Jo 

X J e--cot(0/2) ^ 



JV 

G(0, v?) =#(0, <p) + iJ2 a i [cot(0/2) e~ i(fi ] j ~ l [sin(0 /2)] JV " 1 a(6', 
i=i 

+ !^ r c /6' 1 sm6' 1 [sin(6' 1 /2)]- Ar - 1 [sin(6'/2)] JV - 1 
2 Jo 

v ( 2 \l Wi,<Pi) a(0i,y?i) 

7 ^ 1 cot(0 1 /2) e -«*i - cot(0/2) e"* ^ ibCj 

In the case N = 0, the generalized grand partition function (2.12) reduces to the 
known expression^ for the grand partition function of the two-component plasma. In 
the other limiting case, ( = 0, corresponding to the one-component plasma, the resulting 
expression for the generalized partition function is new ; we show in the Appendix how 
to use (2.12) in this case to reclaim the result of ref. 2 for the exact evaluation of the 
partition function. In the general case N ^ 0, £ 7^ 0, we have not been able to compute 
the grand potential directly from (2.12), but we shall obtain it from a density computed 
in next section. 

3. EVALUATION OF THE CORRELATION FUNCTIONS 

The fully truncated n-particle distribution functions can be obtained by functional 
differentiation of the logarithm of (2.12). Analogous to the situation in the two-component 
plasma limit/ 5 '' the distributions can be expressed in terms of functions which play the 
same role as the Green functions of ref. 5 : 

G S1S2 (0 2 , <p 2 )) := 5 S1S2 5 (($,, (0 2 , <p 2 )) - A S1S2 ((0 l7 Vl ), (0 2 , <p 2 )) (3.1) 



where si,S2 = ±, S((9\, cpi), (02, ^2)) denotes the Dirac delta function on the surface of 
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the sphere 



R 2 Maine dip H(9, (p) S ((0, <p), (0', <p')) = H(9',ip') (3.2) 
Jo Jo 

for any continous function H(9,<p), and A SlS2 is defined in terms of X -1 by 



X" 1 



"ATi a NN F N (9 2 ,ip 2 ) G N (9 2 ,(p2) 

fi(0u<Pi) fN(e lf(pi ) A ++ ((9 1 , lfl ),(9 2 , l p 2 )) A + _((9 1 , lfl ),(9 2 , l p 2 )) 

l9i(0i,<Pi) 9n(0u<Pi) A_ + ((9 u < Pl ),(9 2 ,<p 2 )) A__((9 1 ,< Pl ),(9 2 ,<p 2 )) 



The one-body density of particles of sign s is 



Ps(0,<p) = G aa ((0,<p),(6,<p)) 



(3.3) 



(3.4a) 



the two-body truncated density is 

Pl,s 2 Wi, <Pi), W2, <P*)) = ~G Sl ,s 2 ((01, <pi), (02, <p 2 )) G S2 , S1 ((02, <P2), (0i, <Pi)) (3-46) 

and similarly for the fully truncated n-particle distribution. 

Due to the rotational invariance of the correlations, it suffices to calculate pJ lS2 with 
one of the particles fixed at the south pole (pJ lS2 will then depend on only). Thus, from 
(3.4), it suffices to calculate 



G Sl s 2 ((0, (f), (tt, p>')) := G S1S2 (0, tp) . 



(3.5) 



From the definition (3.3) of the A SlS2 as elements of the inverse of the operator X the 
definition (2.14) of X, and (3.1), it follows that the functions G_ + (0, </?) and G ++ (0,(p) 
satisfy the coupled equations 



/»7T p2lY I 

J d9sm9j ^[sin(0/2)] JV - 1 [cot(0/2) e i,fi ] j - 1 G ++ (6,<p) = I R 2 



7^ J = l 



0, j = 2,..-,N 

(3.6a) 
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G-+(e 1 ,<p 1 )+ % -^ [ d9 2 sme 2 [sm(e 2 /2)] N - 1 (sm9 i y 
^ Jo 

' ^ 2 cot(^i/2) e^i -cot(0 2 /2) 2 # cot(0i/2) ' 



l-AT 



JV 

i 5^ [cot(0i/2) e"^ 1 ]'" 1 [sin(6'i/2)] 7V - 1 G' J (7r) 
i=i 

+ ^ r^ 2 sin^[sin(^/2)]-^- 1 [sin(^ 1 /2)]^- 1 T ^ 2 G "+^ 2 ) 

* Jo Jo 



j Q cot(0 2 /2) e"^ 2 - cot(0i/2) e"^ 1 
+ G++((9 1 ,^i) = (3.6c) 

Also, G_| (0, if) and G (0, (p) satisfy the same coupled equations, with the replacements 

r.h.s. of (3.6a) — > each j = 1, • • • , AT, r.h.s. of (3.6b) — > 
, h ., of (3.6c)^-»l^ e * (3.7) 

Let us consider the coupled equations (3.6) for G |_ and We observe that these 

equations permit a solution of the form 

G ++ (9,<p) = 1++ (9) (3.8a) 

G_+(0, ¥ >)=7-+(0)e-* (3.86) 
G J (n) = , j = 2,---,N (3.8c) 

Substituting (3.8) in (3.6), changing variables t = cos 9 in the resulting equations, and 
setting 

- (1 _ t) t + a + t)1 / 2 (»•*» 
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then gives the coupled equations 

2- (7V " 1)/2 dt 2 h ++ (t 2 ) = (3.10a) 

h-+(t) + 2iTi(R [ dt 2 h++(t 2 ) = 2 (1+N ^ 2 (3.106) 
J -l 2K 

-i 2-( N -^ Gl (n) + toiCRj* dt 2 {1 + + (1 - t)-^h ++ (t) = (3.10c) 

Differentiating (3.10b) and (3.10c) with respect to t, substituting the first of the 
resulting equations in the second, and simplifying gives 

(1-t 2 ) h'L+(t) + (N-l)(l + t) h'_+(t) - (2tv(R) 2 h-+(t) = (3.11) 

This is the Gauss hypergeometric differential equation, in the variable 

s = (l-t)/2 (3.12) 

with (in standard notation) 

a = -l (iV-(jV 2 -4(27rCR) 2 ) 1/2 ) , b = ~ (n + (N 2 - 4(2ttCR) 2 ) 1/2 ) , c=l-N . 

(3.13) 

We note from (3.10c) that we require the solution of (3.11) such that 

h-+(s) = o(s N ~ 1 ) as s^0 (3.14) 

(otherwise the integral in (3.10c) does not exist ; furthermore (3.14) ensures that (3.10a) 
is satisfied). The solution of (3.11) with this property is* 11 ) 

h-+(s) = A s N F(a + N,b + N; N + 1; s) = A s N (1 - s) F(l - a, 1 - b; N + 1; s) (3.15) 
To determine A, we note from (3.10b) the requirement 



/*_+(* = 1) = ^ 2l+f ( 3 - 1G ) 
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This gives 

A _ *C 2l+ ^ r(i - a )r(i - 6) 

"21^ r(iv + i) ^ i7j 

Now that ft |- is known, the value of h ++ follows from differentiating (3.10b) with 

respect to s, which gives 

1 AiVs^ -1 F(a + N,b + N; N; s) (3.18) 



47TZC.R 

Reverting back to the original variable, using (3.12) with t = cos6>, and then to the 
original Green's function, via (3.9) and (3.8), we thus have 

G. + (0, <p) = e-^ r(1 ~^ (1 1} ~ b) cos(fl/2) sin>/2) F(l - a, 1 - 6; N + 1; sin 2 (#/2)) 

(3.19a) 

and, after a simple manipulation on the V functions, 

G ++ (9, <p) = nC 2 r{N+ ^ { ^ + b) sin^- 1 (^/2) F(a + N,b + N;N; sin 2 (0/2)) (3.196) 

To solve the equations (3.6) with the replacements (3.7), we try for a solution of the 

form 

G—(9,ip) = 7— (0) (3.20a) 

G+_(0, ¥ >) = 7 +-(0)e* (3-206) 
F J (n)=0 j = l,---,N (3.20c) 

Proceeding as in the calculation of G ++ and G |_ above, we find 



G+-(e,ip) = G- + (e,ip) (3.21a) 

and 

x F(a + N + l,b + N + 1; N + 2; sin 2 (#/2)) (3.216) 
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4. ALTERNATIVE EVALUATION OF THE CORRELATION FUNCTIONS 

In this approach, we put on the sphere a uniform background charge density of total 
charge — Nq and a point charge Nq at the north pole, with N fixed. Then, we add a variable 
number N + M of pairs of positive and negative particles. We expect that N negative 
particles will stick to the north pole, and therefore there will remain the background, N+M 
positive mobile particles, and M negative mobile particles, as in the previous Section. 

Now, we can use the grand-canonical formalism, with the same variable number N+M 
of positive and negative particles. The charge Nq at the north pole generates a one- 
body potential, which can be taken into account through a charge and position dependent 
fugacity. From (2.1), when V := q 2 /k B T = 2, this fugacity is of the form ([sm 2 (6/2)] sN 
for a particle of sign s(s = ±1). 

The problem on the sphere can be transformed into a problem in the plane by the 
same stereographic projection as in ref. 7. Let P be the stereographic projection of a point 
M of the sphere, from the north pole, onto the plane tangent to the south pole. In terms 
of the spherical coordinates (6>, ip) of M, the complex coordinate of P in the plane is 



The projection is a conformal transformation, with a conformal weight 



That is, an element of length d£ at M and its projection of length \dz\ at P have a ratio 
d£/\dz\ = e w . Also, in terms of the coordinates z and z' of the projections of the particles, 
the pair potential (2.1) becomes 



z = 2R e lv cot - 



2 



(4.1) 




(4.3) 
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Since each particle interacts with N + M — 1 particles of the same sign and N + M 
particles of the opposite sign, the u term in (4.3) gives a total one-body potential q 2 u/2 ; 
for T = 2, this generates a factor e~ w in the fugacity. Furthermore, in the computation of 
the partition function, the area element dS on the sphere can be expressed in terms of the 
corresponding area element d 2 z on the plane as dS = e 2u} d 2 z. Altogether, the system on 
the sphere with a fugacity (e U}sN and the pair potential (2.1) is equivalent to a system in 
the plane with a fugacity 

U\z\) = (e^ sN+ ^ = ( (l + ^) (4-4) 

and a pair potential — qq' In \z — z'\. The N negative particles stuck at the north pole of 
the sphere are projected at infinity on the plane. 

The formalism for dealing with a plane system with a fugacity dependent on the 
charge and the position has been set in ref. 5. The n-particle densities can be expressed in 
terms of Green functions 3 g SlS2 (i"i, r 2 ) (here Vj stand for a position vector in the plane ; 
Zj = rj e llpj ) ; in the plane, the densities are 27r(e' i;(r )(7 ss (r, r), the two-body truncated den- 
sities are — siS2(27r£) 2 e^ ri ) + ^ r2 )|g SljS2 (ri, r 2 )| 2 . They are the stereographic projections 
of densities on the sphere, which therefore are 

p s = (27rC)e-^<7 ss (r,r) (4.5a) 

pL 2 ((0i,^i),(02,¥>2)) = - S i S 2(27rC) 2 e-^)-^)|a SlS2 ( ri ,r 2 )| 2 (4.56) 

where (9j, <pj) is that point on the sphere which is projected at rj. It is enough to fix one 
particle at the south pole and to consider only o SlS2 (r, 0). 

For the present problem, eq. (2.18a) of ref. 5, with now m(r) = 2/rCe w ^ r \ becomes 

(2ttC) 2 e" (r) + A+ e'^A] g ++ (r, 0) = 2yrCo(r) (4.6) 



3 In Section 3 of the present paper and in ref. 5, the Green functions have been defined 
with different irrelevant phase factors, which do not affect the physical quantities. 
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where 



. iuD [ d id N .. , 
or r dip 2 



(4.7a) 



A+ = e~ iv 



_d_ + i_d_ _ N^, 
Or r dip 2 



(4.76) 



Looking for a solution of the form g ++ (r, 0) = o++(r) and changing to the variable s = 
sin 2 | = e w(r ) and to the function k(s) = s~ N / 2 g ++ , one obtains from (4.6) 



d 2 h dh 



(4.8) 



(4.8) is the hypergeometric differential equation. The two solutions are^ 11 -* F((N+5)/2, (N— 
6)/2; N; s) and F((N+5)/2, (N-5)/2; 1; 1-s), where <J = (A^-ie^C 2 !? 2 ) 1 / 2 . The second 
solution behaves like s~ N+1 as s — > ; it would generate in /0+ + (0, 7r) = — (27r(') 2 s Ar ~ 1 |/c(s)| 2 
an unacceptable singularity s _Ar+1 as s = sin 2 (6>/2) — > 0, and it must be discarded. There- 
fore 

g++i r, 0) = C +S ^F^,^;N^ (4.9) 

From (4.6), near the south pole (r = 0, s = 1), o++(r, 0) behaves like — ( lnr ~ — (C/2) ln(l— s). 
The F in (4.9) does have a logarithmic term since 



,'N + d N-5 nT 



T(N) 



r(^)r(^) 



X 



-27 - v> 



N + S 



$ \ ) - ln(l - s) 



as s — > 1 



(4.10) 



where T is the gamma function, •0 is its logarithmic derivative, and 7 is Euler's constant. 
The appropriate behaviour is obtained by choosing 



C+ = 



rgg)rgg) c 
r(iv) 2 



(4.11) 



17 



Then 



<7++M) 



-27 - ip 



N + S 



'(^)-Mi-«) 



as s -> 1 (4.12) 



By a similar calculation, from eq. (2.18b) of ref. 5, one finds 

9 _4,,0) = C_.*+Wl + ^,l + ^;tf + 2;. 



(4.13) 



where 



C_ = 



r(i + ^)r(i + ^) c 
r(A^ + 2) 2 



such that 



g—(r,0) 



-27-^U + ^) - V (l + ^ ) " Ml " «) 



Finally, from eq. (2.18c) of ref. 5, one finds 

N+S\ (N-S 



(4.14) 



, as s — > 1 

(4.15) 



<7- + (r, 0) = ( " + ] ^ ] ^ ni-s 2 ) 1/2 s (N+1)/2 F (l + *±*, 1 + AT + 1; a ) 

(4.16) 



Using the Green functions (4.9), (4.13), (4.16) in (4.5b), one finds 



p++(M) = 



r ^ r ^ . N ^^fN + s n-s Ar . 2 0\ 

ttC 2 V ^tttt sm^" 1 -F , ; iV; sin 2 - 

s T(N) 2 V 2 2 2/ 



(4.17a) 



< 2 



r(i + ^)r(i + ^H) 



t(n + 2) 



sin 



N+l 



xF 1 + 



iV + o , AT-o 



■J + ;iV + 2;sin 2 - 



(4.176) 



C r(i + ^)r(i + V) o. N e w (^ 
2R fWTi) cos 2 sm 2 F (, 1 + 



N + o AT-o 







,1+ — — ;AT + l;sin z - 



(4.17c) 
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Thus, we retrive the correlation functions derived in Section 3 eqs. (3.4b), (3.19), (3.21). 

It can be easily checked that the limit ( — > reproduces the one-component plasma 
results of ref. 2 and the limit N — > reproduces the two-component plasma results of 
ref. 7. 

The thermodynamic limit giving a plane system can also be studied : N, R — > oo, 
6 — > 0, with fixed background density 77 := N/AtiR 2 , fixed fugacity £, and fixed distance 
r = 2Rcos(9/2). By using one of Kummer's relations ( n ) 

z 



F(a, b;c;z) = (1 — z) F c — a, 6; c; 



z - 1 



(4.18) 



iV 



one can rewrite (4.17a) as 

pl + (r) ~ - U r( ^~^ r(a) e-^^iV"( 7 rr ? r 2 )-F f a, a; AT; 1 - 

where a = i^C h- In the limit N -> oo, iV a r(iV - a)/T(N) -> 1 and 



(4.19) 



(yr?7r 2 )- a F ^a, a; AT; 1 - -> e^V^r 2 ) -1 / 2 _ Qj0 (7r?7r 2 ) (4.20) 

where is a Whittaker confluent hypergeometric function/ 11 ) Thus, in the thermody- 
namic limit, 



pl+(r) = - \ V T(a + l)(7r V r 2 )- 1/2 ^i- a ,o(V 



(4.21a) 



Similar calculations give 

p T __{r) = - [vaTia + l^r 2 )- 1 / 2 W_i_ a>0 (^ 2 )] 2 (4.216) 
p^+(r)=p^_= ?7« 1/2 r(a + l)(7r?7r 2 )- 1/2 W_ aj i(7r?7r 2 )] 2 (4.21c) 
The results (4.21a) and (4.21b) had been previously obtained 4 in ref. 10. 



4 The sign of the background charge density has been chosen positive in ref. 10, 
negative in the present paper. Furthermore, there are sign inconsistencies at the bottom 
of p. 127 of ref. 10. 
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5. THERMODYNAMICS 



In the present generalized grand-canonical formalism, the basic thermodynamic func- 
tion is the generalized grand potential Q obtained from the generalized grand partition 
function : 

tt = -k B T\nZ (5.1) 



From (2.4a), the density of negative particles is 

< M > 1 



AnR 2 8nR 2 * d( 



(5.2) 



Thus, we can obtain InS and Q by integrating p-/( with respect to (. 

The densities p s , given by (4.5a), should be independent of r and it is convenient to 
compute them at r = 0. However, for point particles, g ss (0, 0) is a divergent quantity, as 
already noticed in simpler cases^ 5 ' 7 ). For obtaining finite densities, we assume the particles 
to be small hard disks of diameter a rather than point particles. This regularization^ 5 ) 
amounts to replace (7 SS (0,0) by g ss (a,0). Since a is small, we can use (4.12) and (4.15), 
with 1 — s = a 2 /4R 2 , in (4.5a), which gives the densities 



P+ = ttC 



. N + S\ , fN-5\ rtl 2R 
-27-^1— W (— +21n- 



P- = 



-2 7 - V 1 + 



N + 5 



N-5\ rtl 2R 



(5.3a) 
(5.36) 



Using + x) — ip(x) = 1/x, one easily checks that p + — p- is equal to the background 
density rj = N/AuR 2 . 

From now on, we only consider the case of a large system, using large- R expansions 
for fixed values of the fugacity ( and the background density 77 = N/A11R 2 . As R becomes 
large, an expansion of the tj; functions gives 



P- = p + — 1] = rja [—27 — ip(l + a) — Innrja 2 ] + 



a 



AnR 2 



+ a — a ip'(a) 



(5.4) 
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where, again, a = 7r( 2 /r). 

For ( = 0, S should become the partition function of the one-component plasma^ 2 ) 
(the thermal de Broglie wavelength has been taken as unity) 5 

N-l 

Nl 



(5.5) 



p=0 



as confirmed in the Appendix. The corresponding free energy O(0) = —f3 1 lnZ can be 
expanded^ to give 



i i fSJ (0) = 4^|.„^ + iln 



(*™l) 1/2R ] +^-2C'(-1)+o(1) (5.6) 



where (' is the derivative of Riemann's zeta function. 

Using (5.4) in (5.2) and integrating from £ = gives the large- R expansion 



:0 



O(0) + AttR 2 t} [a(27 + Imrqa 2 ) + lnT(l + a)] 



a a 



— — + I x ip'(x)dx + o(l) 



(5.7) 



(kBT) 1 does exhibit the expected universal finite-size correction^ 6 ) (1/3) \nR. That 
correction comes entirely from O(0). 
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5 Our definition of Z includes the usual 1/Nl factor which is not present in the "excess 
partition function" defined in ref. 2. 
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APPENDIX 



When ( = and a = b = 1, 



S 2 (a,6) = Z 2 



(Al) 



where Z 2 denotes the canonical partition function of the one-component plasma on a 
sphere. Using (2.12) we thus have 



Z 2 = An2 Y\ 



(A2) 



where the product is over all eigenvalues of the operator X. The operator X again acts on 
vectors (2.13), except that the component f(6,(p) is no longer present. The eigenvectors 
are therefore of the form 



(k) 
I TV 



and from (2.15) the eigenvalues and eigenvectors are specified by the equations 
Ai-a^ = i d9 sin 



(j = N) and 



JV 



i=i 

We look for solutions of the above equations of the form 



(A3) 



/ dBsva.9 ( dip(sm6/2) N - 1 (cot6/2 e^y^g^ie.ip) (A4) 
Jo Jo 



(A5) 



iktp 



(A6) 



Substituting (A. 6) in (A. 4) gives 

a? = 



(A7) 
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X k a) k) =2nij d0sm8(sme/2) N - 1 (cot6/2y- 1 g k (e) , k=l-j , (A.8) 

and use of these equations gives from (A. 5) 

(A fc - 1) <fc(0) = , fc^{0,---,-(JV-l)} (A9a) 

and 

(l-A fc )A fc = 27T J d9 sin (sin 9 /2) 2(JV_1 )(cot /2)" 2fc 

= 4.lk±j|EW) , * = ,-l,.,-(*-l) . (A96) 

Eq. (A. 9a) gives the eigenvalues = 1 for all k £ {0, • • • , — (N — 1)}, while (A. 9b) gives a 
pair of eigenvalues and A^ for each k = 0, • • • , — (N — 1) such that their product X k \ k 
is equal to the r.h.s. of (A. 9b). 

Thus, writing p = —k, we have 

N-l 

Z 2 = An2 ] I 

p=0 

which is precisely the result of ref. 2, derived from the Vandermonde formalism. 



An 



t{n -p) r(p+i) 
r(N + i) 



(A.10) 
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